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Formulas are derived for the accurate calculation of the Jowest-mode, lumped-element 
representation of perfectly conducting half-round inductive obstacles in rectangular wave- 
guide. These obstacles consist of either one oi two opposed semicircular cylindrical indenta- 
tions extending across the narrow sides of the waveguide. They seem especially suitable for 
use as precise calculable standards of reflection or impedance in waveguide. Schwinger's 
integral equation approach [1,2] 1 is used to obtain stationary expressions for the desired 
parameters as functionals of the surface currents on the obstacles. Upper bounds are ob- 
tained for one of the two parameters. Explicit formulas are derived for the values of the parameters 
under the assumption of n-term Fourier sine-series expansions for the obstacle currents. Rapid 
convergence is indicated by numerical evaluations for n=l, 2, and 3. In the process of 
obtaining expressions suitable for numerical calculation, an expansion (believed to be new) 
of the Green's function of the problem is obtained and the sums of certain infinite series of 
Bessel's functions occurring in this expansion are expressed in terms of definite integrals. 
A brief numerical table of these sums, sufficient for the evaluation of the n = l approximation, 
is included. 

1. Introduction 

In this paper formulas are derived for the accurate calculation of lowest-mode, lumped- 
clement parameters for what may be identified as "single half-round" and "double half-round" 
inductive obstacles in rectangular waveguide. As shown in figure 1, the obstacles consist of 
semicircular cylindrical indentations extending across the narrow sides of the waveguides. 

The particular geometry considered seems especially suitable for obstacles to be used as 
calculable standards of reflection or impedance in waveguide. The geometry is well suited to 
clectroforming, so that obstacles may be fabricated by this process as well as by machining. 
For obstacles producing standing-wave ratios of moderate values, the obstacle radius is large 
compared to high-standard machining tolerances, so that unduly close tolerances are not 
required. A change in waveguide cross section before and after the obstacle is avoided, elim- 
inating the need for waveguide components in odd sizes and permitting maximum flexibility 
of interconnection. 

In the present work the obstacle and waveguide surfaces are assumed perfectly conducting. 

The Waveguide Handbook [3] contains approximate formulas for semi-elliptical obstacles, 2 
which can be specialized to apply to the present problem. However, for the contemplated use 
in standards work, formulas with very much greater accuracy (10 3 or 10 4 times greater) are 
wanted. 

Brief descriptions of the present and some related work have been given previously [4]. 

2. Formulation of Problem 



We employ the customary complex electric and magnetic field vectors E=E(r), H=H(r), 
which satisfy Maxwell's equations in the form 

VXH=jueE+J, [ 



1 Figures in brackets indicate the literature references at the end of this paper. 

2 The pertinent formulas are afflicted with a number of misprints; correct forms of these results for half-round obstacles are given later in this 
paper. 
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under the assumption of time-dependence represented by the (omitted) factor exp ( jcd) . Here 
t is the time, u/(2w) is the frequency, j is the imaginary unit, / is the complex current-density 
vector, and /z and e are respectively the permeability and the permittivity of the homogeneous, 
isotropic, nondissipative medium in the waveguide. The MKS system of units is employed. 

As shown in figure 1, we choose a rectangular coordinate system Oxyz such that the in- 
terior of the waveguide is the space 0<£<a, 0<y<J, — oo<£<oo and the surfaces of the 
obstacles are given by the loci 

R 2 =x 2 +z 2 , 0<x<R, 

- - y (2.2) 

R 2 =(a-x) 2 +z 2 , 0<a-x<Rj 

in the double half-round case and by 

R'=x 2 +z 2 , 0<x<R (2.3) 

in the'single half-round case. 



JO 





JO 




Figure 1. 



Single and double half-round inductive obstacles 
in rectangular waveguide. 



The fields incident on the obstacles are, by hypothesis, to be in the TE 10 mode. Traveling- 
wave forms for this mode may be written 

E±=sin (Tx/a)e^e y , 1 

H ± = [Trj sin (TTx/a)e x +JTr(ccfjia)~ l cos (wx/a)e 2 ]e Tj P z ' j 

The upper and lower signs refer respectively to waves going in the positive and the negative 
^-directions; the phase constant /3 is equal to [k 2 — (ir/a) 2 ] l/2 , where & = co(jue) 1/2 ; the wave-ad- 
mittance 77=0/(co/x); and e x , e y , e z are the unit vectors of the s^ystem Oxyz. It is assumed that 
7r<te<27r, so that the TE 10 mode, but no higher mode of the type TE re0 , will have real propaga- 
tion. 

We define 'Voltage" and "current" v m {z m ) ,i m (z m ) for the TE i0 mode by means of the 
equations 

E my =v m (z m ) sin (irx/a), I 

fm=l,2 (2.6) 

H mx =( — ) m i m (z m )r) sin (irx/a).j 

Here the indices m=l,2 refer respectively to the "left-hand" side (zi< — B) and the "right- 
hand" (z 2 y > R) of the obstacle considered, and E mv , H mx are components (identified by the 
subscripts) of the TE 10 -mode part of whatever total electromagnetic field may be present in 
the waveguide at z = z m . From (2.6) and (2.4) it follows that the waveguide characteristic 
impedance, defined as the value of v m (z m )/i m (z m ) for a pure traveling wave incident on the side 
m of the obstacle, is equal to unity for m = l,2. 



114 



The desired lowest-mode, lumped-element description of an obstacle is contained in its 
impedance matrix Z, which characterizes the linear relations 



Vi = Znt'i + Z 12 z 2 , » 2 = Z 2 it'i + Z 22 i 2 



(2.7) 



imposed by the obstacle on the terminal variables v m ,i m . The values of the Z ron depend upon 
the location of the terminal planes z=z x and 2=£ 2 ; in the presenl problem it turns out to be 
convenient to extrapolate these planes to 2i=z 2 =0. 

Now, the equality Z n =Z 2 2 follows from the structural symmetry (and the symmetrical 
disposition of the terminal planes) with respect to the plane 2=0; Z 12 = Z 2 i is assured by both 
symmetry and reciprocity; and, since losslessnoss has been assumed, Z is pure imaginary. 
Thus Z is of the form 



ljX n ]X a ,\ 

1./-V,, jX n )' 



(2.8) 



and there are only two independent parameters to bo determined. 

We do not obtain formulas for X n and X l2 directly; but rather, to exploit the symmetry 
of the problem more fully, we consider the impedances obtaining under modes of excitation in 
which the electromagnetic field is either symmetric or antisymmetric with respect to the plane 
2=0. Thus in the symmetric case, p 1 (0)=t? 2 (0), &i(0) = / L >(0). and the impedance may be 
written 

Z-=f>i,.(0)/ti ># (0), (2.9a) 

where we have added a subscript distinguishing the "even" case. In the antisymmetric 
("odd") ease, /*,(()) — v 2 (0), &i(0) — 2 2 (0), and the impedance may be written 



Zoo=ci,o(0)/ii,o(0), 
again adding a distinguishing subscript. The relations 

2jX, , = Z ee + Z ()0 , 2jX l2 = Z ee - Z 0( 



(2.9b) 



(2.10) 



follow directly from (2.7), (2.8), and (2.9). (It may be noted that the symmetric and anti- 
symmetric field distributions correspond to the eigenvectors of the matrix (2.8) and the quanti- 
ties Z ee and Z 00 are the eigenvalues of this matrix.) 

An equivalent network for the obstacles is of interest and is presented here for convenient 
reference. For a V-notwork, using the si<m conventions shown in figure 2. we find 



Z\ — Z(ki. 2Z 2 — /j ee Z o, 



(2.11) 



where Z, and Z 2 are respectively the series and the shunt elements of the symmetrica] T. It 

may be remarked that the obstacles considered are called "inductive" because for small radii 
tin* important element, Z 2 , is a positive reactance. 




Z,= It 



O (Zap ~ Zaq) 



' *-00 L 2~z Kt -tt~ *-00' 

Figure 2. Equivalent network. 



54253S <',<> 
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3. Analysis 
3.1. Reduction to Scalar Form 
Maxwell's equations directly imply the general differential equation for E, 

VXVXE=k 2 E-jawJ, 



(3.1) 



which must be specialized appropriately for the present problem. E is moreover subject to 
the boundary condition that its tangential components vanish on perfectly conducting surfaces. 
In the present problem both the exciting field and the waveguide-obstacle structure are 
independent of y in the range 0<?/<&, and it follows that the whole field, secondary as well as 
primary, will be independent of y in the same range. Further, since the x- and .^-components 
of E must vanish in particular on the y = surface of the waveguide, these components must 
vanish for all values of y involved. Thus the electric field is of the form E=<j>{*\ z)e y ; eq (3.1 ) 
(with /=0) reduces to 

V 2 <£+& 2 <£=0 in S, 

and the boundary condition becomes 

</>=() on (\ 

where S is the cross section and C is the boundary of the structure in a plane y= constant (fig. 3). 

x 




Q S /( 



Figure 3. Relevant geometry and coordinates for the single 
and double half-round problems. 

The magnetic Held corresponding to E of the above form is 

.-<«-(g*-g«,) 

The surface current K=HXn on C is thus given by 



(3.2a) 



(3.2b) 



evaluated on C, where n is the outward normal unit vector on O. Thus K has only a //-com- 
ponent and may be written as K(x,z)e y . 

Inasmuch as the whole problem can be dealt with in terms of c/> and K, the problem is es- 
sentially a two-dimensional scalar one. In what follows it will be convenient to refer to 
simply as "the field". 
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3.2. Integral Equation Formulation 

We now introduce, a Green's function T=T(x, z, x' ', 2') such that r is the field that would 
be produced in the unperturbed waveguide at the point (,r, z) by a unit current filament par- 
allel to e v at the point (V, 2'). As a function of (x, z), T satisfies the differential equation 

v 2 r+Fr=ico M 8(x-x') s{z-z f ) (3.3) 

obtained from (3.1) with J=b{x — x r ) b(z—z f )e v , where <5 denotes the Dirac delta-function. 
T further satisfies the boundary condition r = for x = 0, a and represents outgoing waves 
for \z— 2'|— >oo. An explicit form for r is [5] 

T(x, z, x',z')=—a)fx{pa)~ 1 sin (irx/a) sin (wx'/a)e-M z -^ 

— (joon/a)^,^ 1 sin (nirx/a) sin (n,Tx'/a)e~ anlz ~ z '\ (3.4) 

n = 2 

where a n =-\ (nw/a) 2 —k 2 for n>2 and the other symbols are as already defined. 

It is convenient to put 

T=-MF+jG), (3.5) 

where F and are real. We observe that /'can be written in the form 

F=(Pa)- 1 lfe(x, z) *.(*', a') +*>(*, 2) *o(*', *')], (3.6a) 

with the definitions 

^ e (a, 2) = cos /32 



^ (», s)=sin ^ 



te sin(ira?/a), 
,2 sin(7rj:/a). 



(3.6b) 



The functions ^ e and ^ are respectively symmetric and antisymmetric with respect to re- 
flection in the plane 2=0. 

With the aid of the Green's function we can set up the fundamental equations of the prob- 
lem in integral-equation form. Let \p=\j/(x, z) denote the exciting field i.e., the field, due to 
remote sources, that would exist in the absence of obstacles. Then, with an obstacle present. 
the total field </> can be written as the sum of the exciting field and the secondary field due to 
the current on the obstacle : 

4>(x,z)=1>(z,z)+ f T{x,z,z* ',z')K(x' >') dC, (3.7) 

where K is the surface current on the obstacle and the integral is a line integral going over the 
contour C of the obstacle (C consists of two parts in the case of the double half-rounds). If 
the point (a?, z) is taken on Co, </> must vanish and we have 

0=f(x,z)+( T(x,z,x',z')K(x',z')dC, (x,z) on C . (3.8) 

This is an inhomogeneous integral equation of the first kind determining i£for a given \p. The 
equation reduces to slightly different forms in the symmetric and antisymmetric cases, which 
we now consider. 

a. Symmetric Case 

Let the exciting field be A\l/ e (x, z), where A is an arbitrary amplitude and \// e , as defined in 
(3.6), is a symmetric function of z. The surface current on the obstacle will then also be a 
symmetric function of z; we denote it by K e . We now examine (3.7) for 2<<V, say, to find 
the lowest-mode component and thus to find Vi, e (z) and i Ue (z). For 2<<V, (3.7) becomes 

<t>(x,z)=A cos pz sin (wx/a)—(pa)- l <ap sin (wx/a)e m I ^ e K e dC, (3.9) 



.vh.t.::n <»<) 
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since integrals of the type \poK e dC vanish because of the mutual orthogonality of sym- 

J c o 
metric and antisymmetric functions. By comparison with (2.6) it is seen that 

v he (z)=A cos fiz-{fia)~W e^ z f $ e K e dC. (3.10) 

J c o 

Using (3.2a) to obtain the ^-component of the magnetic field associated with (3.9) and again 
referring to (2.6) one finds 

iiM)=-JA sin jte+Gto)- V e* 1 f f e K e dC (3.11a) 

Equations (3.10) and (3.11a) hold directly from the definitions for z<^ — R and by extrapolation 
for z=0. In particular (3.11a) yields 

ii..(0) = (jte)" 1 «n (teKedC. (3.11b) 

Next, in the integral eq (3.8), we separate the Green's function into real and imaginary parts as 
in (3.5), use the orthogonality property of even and odd functions, and thus find 

fc(x,s) v ltB (p)=Jw \g{x,z,x'z>) K e {x',z f ) dC. (3.12) 

This equation and (3.11b) together furnish a definitive mathematical statement, of the sym- 
metric part of our problem. 

b. Antisymmetric Case 

In this case we let the exciting field be A\f/ (x, z), where A is again an arbitrary constant and 
^o, defined in (3.6), is an antisymmetric function of z. The current, K , on the obstacle will 
then also be antisymmetric with respect to z. In the same manner as in the symmetric case 
one finds, for z<0, 



In particular, 
Further, 



Vl (z)=A sin /3z + (/3a)-y« M g# 2 f toK dC. (3.14a) 

»i.o(0) = (j8a)- 1 i«/*f foK<4C. (3.14b) 

Jc 



ix *(z)=jA cos 0z — (Pa)- 1 jane** \ toK^C. (3.15) 

and the integral equation becomes 

*»(*,2)<i.o(0)=-a V i f G(x,z,x f ,z')Ko(x',z')dC. (3.16) 

Jc Q 

This equation and (3.14b) together furnish a definitive statement of the antisymmetric cases 
of our problem. 

3.3. Summary of all Cases 

In the case of the double half-rounds the above integrations go over a two-part contour 
consisting of an upper semicircle C u and a lower semicircle C h An equivalent problem involving 
only one of the parts of the contours, say C h can be formulated using the fact that the exciting 
field as well as the structure is symmetric with respect to the plane x=a/2, so that the current 
on the obstacles must also be symmetric with respect to this plane. This symmetry is quite 
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independent of the symmetries with respect to z=0 distinguished by the subscripts 'V' and 
"o," and so these subscripts will be dropped for (he moment. Thus, it is clear that 

( tKJC=( xfsKdC+i xpKdC=2 f ypKdC. (3.18) 

f G(x,z,x',z f )K(x',z')dC'=( G {2 \x,z,x' ,z') K(x' ,z') dC , (3.19) 

J C o J C l 



Further, 



where G (2) is defined by 

G {2 \x,z,x' ,z')=G{x,z,x' tZ^+Gix^a-x' ,z'). 

Equations (3.18) and (3.19) enable the desired restatement of the double half-round problems. 
In what follows, superscripts 1 ,2 will be used when it is desired to distinguish quantities 
associated with the single and the double half-round problems, respectively (this means in 
particular G=G il) ). Integrals will be indicated by means of the convenient scalar product 
notation, e.g., 

ty,K)=f c fKdC, (3.20a) 

and by the operator notation 

GK= f G(x,z,x',z') K{x\z f ) dC. (3.20b) 

J c i 

As indicated, all such integrals are to be taken over C h the lower semicircle. As a further 
notational convenience, we introduce the normalized" sin-face currents, 

L=]»liKJviJP), I o =-wK o /ii,o(0), (3.21) 

for both the single and the double problems. 

Our results thus far may now be summarized as follows. For the symmetric cases we have 

+ e =G<»I e , l/Xl h ?=h(+ e ,I e )K(ki) (3.22) 

and for the antisymmetric cases 

+ = <?<»/„ Xf =-h(+ ,I )/(Pa), (3.23) 

where h=\ or 2 for the single or the double half-rounds, respectively, and X ee , X 00 are the 
reactances corresponding to the (pure imaginary) impedances Z ee , Z 00 defined in (2.9). 

4. Solution for the Reactances 

4.1. Application of Rayleigh-Ritz Method 

In each of the problems specified by (3.22) and (3.23), an unknown reactance or suscept- 
ance is proportional to a scalar product 

M=(+,I), (4.1a) 

where the function / is determined by an integral equation, 

GI=x(s, (4.1b) 

in which \p and G are given and G is symmetric. From these equations one may easily con- 
struct the "stationary representation" [1,2] 

M=(f,I)V(I,GI) (4.2) 
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for M. This constitutes the starting point for the method of Rayleigh and Kitz [6], which 
we wish to use to obtain an approximate expression for M. One assumes in (4.2) an approxi- 
mation for the surface current in the form 

where the x t are coefficients to be determined, the j t are members of a suitable set of basis 
functions (to be chosen explicitly), and summation from 1 to n over repeated indices i,j, . . ., 
is understood. Expression (4.2) becomes an ordinary function of the x t and is to be subjected 
to the conditions of stationarity dM/dXi=0, i=l, 2, . . ., n. This leads to the system of 
equations 

G iJ x j =(x k c fc /M n )c i 
where 

d={+Ji), 0»-(f t ,Qf,), (4.4) 

and M n denotes the now-determinate approximate value for M. If the aibitrary normalization 
of the x t is chosen so that x k c k =M n , then 

Xi=(G- jijCj, 

where {G~ l ) tj denotes an element of the inverse of the nXn matrix of the G ih and it immediately 
follows that 

Mn^iG-^ijdCj. 

Finally, we write this quadratic form as a ratio of determinants, 



-Mn 





Ci 



Ci 

Gn 



Gin 



G n i 



G n 



G n 

G„i 



Gin 



(4.5) 



Thus M n is expressed in terms of the known quantities (4.4) ; if the objective is only to calculate 
the reactances, it is not necessary to calculate the x t . 

4.2. Definiteness of G; Upper Bounds for X 00 

Although in our problem G is not a definite kernel, it does have definite character with 
respect to a suitably restricted class of admissible functions. In fact, G is positive definite 
with respect to functions that are antisymmetric with respect to z and vanish for |^|>X g /4, 
where \ g =2w/l3 is the "guide wavelength" of the single propagating mode. 3 In terms of the 
integrals of interest here, this means that if 7 is any assumed surface current distribution, 
continuous, not identically equal to zero, and antisymmetric with respect to z, then 

(7 , GI ) >0, (4.6) 

provided E<X g /4. 

Now, for a positive definite kernel it is well known that the right-hand side of (4.2) is not 
only stationary, but also a maximum for the true /. For this result we offer the following 
brief proof based directly on the Schwarz inequality. 4 Assuming G to be positive definite, 
we wish to compare the approximate value for M given by 

M={I,WI{I,6I), 



3 This result is suggested by qualitative physical considerations and has been verified analytically by the author. It has also been verified 
that G is not definite with respect to functions symmetric with respect to z. 

* Proofs of this result in the literature known to the author, e.g., [1,2], employ the so-called bilinear expansion for the kernel involved and do 
not yield the "only if" part of the result (4.7). 
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whore I is an admissible approximation for the surface current, with the true value given by 
(4.1a) (as well as by (4.2)) when I satisfies (4.1b). From (4.1b) it follows that (7, ^) = (7, 01); 
the Schwarz inequality for positive-definite symmetric transformations [7] gives 

(7, GI)*<(I,GI) (J, GI). 

Since (/, GI) = (I^)=M, we have the result 

M<M. (4.7) 

Moreover the sign of equality holds in the Schwarz inequality and hence in this result if and 
only if I is proportional to /. 

Clearly, (4.7) is applicable when G is merely positive definite with respect to all functions 
that need be admitted; by (4.6) this will be the case in the antisymmetric parts of our problem 
provided i?<X g /4. Hence, subject to this inequality; in the antisymmetric cases we shall have 
M n <CM (we do not anticipate an exact result for finite n); this in turn implies that numerical 
values for X 00 will be (algebraic) upper bounds for the true values- — assuming, of course, that 
numerical evaluations involved are sufficiently precise. 

Unfortunately one can not obtain a bound for X ee in this simple manner. Furthermore, 
derived quantities of practical interest (such as (4.17)) will usually depend upon both X ee and 
X 00 . Thus in general no bounds for such derived quantities are determined. 

4.3. Basis Functions 
We introduce plane polar coordinates r,6 such that 

X—-T sin 6, z=r cos 6. (4.8) 

The path of integration C t is then given parametrically by 

x=R sin d, z=R cos 0, (4.9) 

where R is the obstacle radius and 0<#<7r. On this path the surface current I=I(R sin 0, 
R cos 0) becomes a function of 0, which we denote simply by 1(d). We observe that sym- 
metries with respect to z=0 are equivalent to symmetries with respect to d=w/2. 
As basis functions we take 

j t = 2(jR)- 1 sin (2i-l)6 (4.10) 

for the symmetric cases (expansion of I s ), and 

j i =2(ivR)- 1 sin 2id (4.11) 

for the antisymmetric cases (expansion of I ). The two sets of functions are complete for the 
expansion of symmetric and antisymmetric functions, respectively, in the interval (0,7r). The 
completeness gives good assurance that for sufficiently large n, M n will approach arbitrarily 
close to the true M. For practical numerical calculation it is vital that the convergence of M n 
be rapid. In the present instance rapid convergence might be anticipated on the grounds that 
the true I e and I may be expected to be smooth (possibly infinitely smooth) functions of 
vanishing at and t, so that their Fourier sine series should converge rapidly. 5 



5 The points (R, 0) and (R, it) are the vertices of the corners formed where the half-round centered at meets the waveguide wall. The van- 
ishing of the surface current in these vertices is equivalent, by (3.2), to the vanishing of the components of the surface magnetic field normal to 
the vertices. That these field components should vanish in the limit as a vertex is approached follows from properties of Sommerfeld's solution 
[14] of the problem of diffraction by a wedge of arbitrary angle. 
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4.4. Evaluation of Integrals 

Integrals of the type i}p,j x ) may be obtained relatively easily with the aid of the well-known 
expansion [8] 

n = — a) 

wherein J n is the n th order Bessel function of the first kind. In this expansion we put £=kR 
and a=±d — x+(W 2 )> with k sin x=V&> h cos x = 0, and employ (4.9). The expansion 



e ipz sin (7rx/a) = 2^ l j n ~ 1 J n (kR) sin nd sin ??x 



(4.12) 



is then easily obtained. The integrations yielding the desired quantities are now elementary; 
one obtains 



c t =2(-y p - 1)/2 J p {kR) sin p x , p=2i-l, 
for the symmetric cases, and 

c i =2{-)^- 2) ' 2 J p (kR) sin p X , p = 2i, 



(4.13a) 



(4.13b) 



for the antisymmetric cases. (It is encouraging that these quantities decrease rapidly with 
increasing index.) 

The evaluation of the integrals of the type (/*, 6/^ = 6^ is considerably more complicated 
than the evaluation of the c i} and the task has been relegated to the appendix The results 
obtained may be summarized conveniently in the following form, which comprehends four 
cases : 



Q(h) 



1 



CiCj 4to p-S q 



[p P (kR)8 pq +a p h i(ka)] 



(4.14) 



Here s^sin^x; p p (kR) = Y p (kR)/J p (kR); in the symmetric cases, p = 2i—l and g=2/-— 1; in 
the antisymmetric cases, p=2i and q_ = 2j; the functions a p h q are defined in the appendix (and 
tabulated briefly in table 2); Y p denotes the Bessel function of the second kind; and all other 
symbols involved have been defined previously. 



4.5. Results 

The value of a reactance element corresponding to an Ti-term expansion (4.3) will be 
called an "n th approximation" and will be written simply XfJ or X$ without special notation 
indicating the value of n being considered. An explicit expression, in determinantal form, for 
Xfe in the n th approximation may be obtained by combining (3.22), (4.5), and (4.14); viz, 



X ih) -- 





Px+vl? . • 


„(h) 
°l,2n-l 




$a 




°l,2w-l • • 


P2n-\~T (T 2n-l,2n-l 




4h 





»1 

Pl + *i? 


• • • s 2n-l 
„(h) 






S 2n- 


1 Ol,2re-l 


• • • P2n -1T%-1, 


m-i 



(4.15) 
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(Iii this and in the following equation the arguments kR and ka are understood.) 

(3.23), (4.5), and (4.14) yield 



Similarly, 



P2+*$ 



°2. 2n 



a 2,2n 



pa P2-TV22 



P2n~y~02n,2n 



°2, 2n 



a 2, 2n 



P2n~\°~2n, 2n 



(4.16) 



Hie determinantal form of these expressions seems convenient enough for practical calcula- 
tions with n at least as large as 3. However, it may be remarked that, because of the special 
form of the diagonal elements, the determinants can be expanded in the same way as certain 
determinants occurring in the Fredholm theory of integral equations. This expansion may 
be convenient for the consideration of arbitrarily large n, but is not needed here. 

Once the reactances X ee and X 00 are determined, other quantities or parameters associated 
with the representation of an obstacle are of course also determined. In particular, the "volt- 
age" standing-wave ratio seen on one side of an obstacle when the other side is terminated 
in a matched (reflectionless) load is a familiar quantity having immediate physical significance, 
and it is convenient to discuss some of the results in terms of this quantity. Recalling that 
waveguide characteristic impedances were chosen equal to unity and referring to the equivalent 
network, figure 2, it is easily found that the standing-wave ratio in question is given by 



„=(1+|S|)/(1-|£|), 



where \S\ , the magnitude of the associated reflection coefficient, is 



(4.17a) 



(4.17b) 



In order to get some idea of the behavior of the sequence of approximations, the reactance 
elements and the related VSWK (4.17) have been calculated for n=l f 2, and 3 and for a range 
of values of ka and kR. Table 1 presents more or less typical results in terms of the values 
of VSWR. The convergence of the sequence of approximations, judged on the basis of 
numerical results exemplified in the table, appears to be very rapid; by the same token, the 
first approximation furnishes a rather good result for VSWR's up to about 2. 

Thus far it has not proved feasible to determine the manner of convergence of the sequence 
of approximations by theoretical means. In view of the character of the particular problem 
at hand, it seems probable that the numerical evidence may be relied upon. 

Table 1. Convergence of VSWR 

ka=4.5 



kR 


Approxi- 
mation 


VSWR, 

Single half- 
round 


VSWR, 

Double half- 
round 


0.2 


1 I 


1.0370968 
1.0870970 

1. 0370970 


1. 0776492 
1.0770499 
1. 0776499 


0.7 


( I 


1.4547608 
1.4554655 
1.4554655 


2.8411681 

2.8410260 
2. 8416268 


1.0 


1 1 


2. 1070625 
2. 1125072 
2.1125112 


15.900756 
16. 009131 
16. 009479 



Note: Maximum computational inaccuracy in 
the above figures is estimated to be approximately 
±1 in the sixth significant figure; eight figures, as 
given, may be significant in indicating the behavior 
of successive approximations. 
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Further calculations have been made regarding ka and R/a as independent variables. 
This choice of variables corresponds to the practical situation where one has a given obstacle 
(having a given value of R/a) that is to be operated at various frequencies. Figures 4 and 5 
show the matched-termination VSWR (4.17) resulting for various values of R/a and f/f c , 
where / is the operating frequency and j c is the cut-off frequency for the TE i0 mode (and 
jlj c =kalir). (Precise tables (based on the third approximations) giving basic parameters and 
possibly some derived quantities are being prepared for separate publication.) 

The first approximations appear to be of useful accuracy for some purposes and represent 
an appreciable extension of previously available results. Consequently, a brief tabulation 
of the o--f unctions, sufficient for the calculation of first approximations, is presented as table 2. 
For convenience of reference we write down (4.15) and (4.16) for the special case n=l } 
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Figure 4. VSWR versus f/f c for single half-rounds. 
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Figure 5. VSWR versus f If c for double half-rounds. 



x™=- 



13a 
4hsl 



[p 1 (kR) + ai h l \ka)l 



X ih) — 

■**. no — 



ihsl 

Pa 



[f>2QcR)+c%\ka)]-\ 



(4.18) 
(4.19) 



and note that Si=w/(ka) and s 2 =%n-p/(k 2 a). 





Table 


2. Values of 


*£(*) 




X 


•fffr) 


* { £(x) 


»i?<w 


»g>fr) 


3.80 
3.85 
3.90 
3.95 


2. 055435 
1.976710 
1.899606 
1. 823995 


-0.3170852 
-.1982432 
-.09078960 
006197600 


3. 960997 
3. 857142 
3. 756554 
3. 659060 


-2. 544410 
-2. 215670 
-1.912736 
-1.633535 


4.00 
4.05 
4.10 
4.15 


1. 749764 
1. 676799 
1.604988 
1. 534222 


. 09356160 
. 172C644 
. 2424036 
. 3052036 


3. 564510 
3. 472757 
3. 383659 
3. 297087 


-1.376168 
-1. 138910 
-0. 9201812 
-. 7185580 


4.20 
4.25 
4.30 
4.35 


1.464400 
1. 395414 
1.327159 
1. 259533 


. 3610476 
.4104564 
. 4539056 
. 4918312 


3. 212923 
3. 131048 
3. 051353 
2. 973735 


-.5327088 
-. 3614372 
-.2036444 
-.05831720 


4.40 
4.45 
4.50 
4.55 


1. 192430 
1. 125746 
1. 059372 
0. 9931968 


. 5246272 
. 5526568 
. 5762488 
. 5956964 


2. 898094 
2. 824341 
2. 752387 
2. 682146 


. 07546600 
. 1985600 
.3117472 
. 4157360 


4.60 
4.65 
4.70 

4. 75 


.9271072 
. 8609836 
. 7946988 
. 7281208 


.6112796 
. 6232416 
. 6318024 
.6371680 


2. 613542 
2. 546500 
2. 480946 
2. 416812 


.5112028 
. 5987500 
. 6789328 
. 7522756 


4.80 
4.85 
4.90 
4.95 


. 6611076 
. 5935064 
. 5251500 
. 4558572 


. 6395192 
. 6390196 
. 6358108 
. 6300236 


2. 354034 
2. 292551 
2. 232301 
2. 173228 


. 8192564 
. 8803180 

. 9358596 
. 9862656 


5.00 
5. 05 
5.10 
5.15 


. 3854276 
.3136376 
. 2402376 
. 1649416 


. 6217692 
.6111440 
. 5982284 
. 5830836 


2. 115279 
2. 058400 
2. 002545 
1. 947661 


1. 031884 
1.073038 
1.110027 
1.143125 


5.20 
5. 25 
5.30 
5.35 


. 08742720 
. 007320000 
-. 07581400 
-. 1624856 


. 5657628 
. 5462944 
. 5246972 
. 5009664 


1.893705 
1.840631 
1. 788397 
1. 736963 


1.172594 
1.198667 
1.221571 
1.241508 


5.40 
545 
5.50 
5.55 


-. 2533004 

-. 3489884 
-. 4504328 
-. 4587224 


. 4750804 
. 4469912 
. 4166292 
. 3838912 


1. 686288 
1. 636332 
1. 587059 
1. 538434 


1.258672 
1.273236 
1.285370 
1. 295225 


5.60 
5.65 
5.70 

5.75 


-. 6752200 
-. 8016548 
-. 9402632 
-1.094004 


. 3486376 
.3106864 
. 2697956 

. 2256488 


1. 490420 
1. 442981 
1. 396088 
1.349701 


1.302940 
1.308654 
1.312488 
1. 314555 


5.80 


-1. 266887 


. 1778340 


1. 303798 


1.314964 



Note: Accuracy of the above figures is estimated to be approximately one 
part per million. 



124 



For sufficiently small kR ) the expressions just given may be simplified still further. For 
small kR, the quantities Pi(kR) and p 2 {kli) become large, in accordance with the formula 



Pt (kR)=*-p\ (2-1)! (j^)* (2^1), 



and thus determine the qualitative behavior of X„ and A r 00 . Using this approximation for 
the p's, neglecting the <r's, and introducing the guide wavelength X g =2ir/j3, one obtains 



Xh 



w = 2o/_ffl_V f42Q v 

!e AX, w?/ 



«=4:(f)' <«•> 

(These are the correct forms of the results in the order of approximation of the more general 
formulas in the Waveguide Handbook [3].) Although these formulas give a rather good picture 
of the behavior of the reactances, they fail to reflect the interesting qualitative feature that 
X™ and X% must eventually become equal to each other as R approaches a/h (as can be seen 
physically). 

5. Appendix 

5.1. Expansion of Green's Function in Terms of Cylindrical Wave Functions 

In order to arrive conveniently at the respective types of Green's functions needed for 
the single and the double half-round problems, we consider a waveguide of width w (instead 
of a), bounded by the lines # = and x = w; and, in addition to the electric boundary condition 
at x = 0, we consider both electric and magnetic boundary conditions at x=w (by "electric" 
or "magnetic" boundary conditions it is meant that the Green's function or its normal deriva- 
tive, respectively, vanishes on the boundary). The radius vectors of the source-point and the 
field-point will be denoted by r' and r, respectively. 

Now, it is well known that tin 1 two-dimensional "free-space" Green's function, satisfying 
(3.3) and representing outgoing waves at infinity, is —(w/x/4)// (/"|r / —r|), where H denotes 
the Hankel function of the second kind (as is appropriate for exp (j<d) time dependence). 
It is also well known that the desired Green's function, <l>, satisfying (3.3) and the required 
boundary conditions at x=0, x=W, and at 2=±°°, can be obtained in the form of a sum of 
free-space Green's functions by the method of images [9]. In fact 

$=-(«m/4)S s"[Ih{k\pt-r\)-Ih{k\p--r\)}, (5.1) 

n = — oo 

where p^=(2nw±x/)e x J rz'e z and s= + l or —1 according to whether the boundary condition 
at x=w is electric or magnetic. By inspection of the array of sources (fig. 6) it can be seen 
that this function has the proper symmetry to satisfy the boundary conditions: it is anti- 
symmetric with respect to x = and it is antisymmetric or symmetric with respect to x = w 
according to whether s=-f-l or — 1. 

Equation (5.1) can be transformed into a sum of elementary wave-functions of polar 
coordinates by repeated application of Graf's addition theorem for Bessel functions [10]. 
We introduce polar coordinates (r, 6), (/, Q f ), and (p tt , B n ) for r, r' , and p n in the manner shown in 
(4.8), and we assume r and r' both to be less than w, so that p n ^>w for n ^ 0. Considering 
the triangles with sides p n , r, and p n — r and applying the addition theorem, we obtain 

4W * P ^„J P (kr')H p (kr) e {e ] 

ft = —oo p= — oo 
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-2w 



Figuee 6. Source points in the method of images. 

Relative signs of sources are indicated by the adjacent symbols +, 
+s, and — s; the role of s is defined in connection with (5.1). 



where the terms arising from n=0 appear in the first line and are to be omitted from the primed 
summation in the second line; the upper and lower alternatives in the first line apply according 
as r<jr' or r^>r', respectively. Again, by the addition theorem, 

H p {k P 1i)e-J<= S J g (kr')H p+g (\2nkw\) exp [j(tt/2) (q-p) sgnn±jq6']- 
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Combining the last two equations and rearranging terms, one finds that (5.1) can l>c written 
in the form 



-t>= — COjlt 



ggsm^smg* | j p(kr >)H,(kr) 8 "« 

+2J p (kr)J,(kr') cos [(t/2)(/j- 2 )] f} «^ff,_ f (2n*tiO-.H, + ,(2nifcw)] "V ■ (5.2) 

71=1 J 

If now one puts w=a and s= + l, this expression becomes the desired representation of V of 
the text; similarly, if w=a/2 and s= — l, (5.2) becomes the desired representation of a Green's 
function r (2) related to 6? (2) in the same way as T is related to G (cf. 3.5). Thus, employing the 
relation (3.5), the relation H p —J p —jY p , the notation (3.20), the basis functions (4.10, 4.11), 
and carrying out the integrations with respect to 6, 6', one obtains 

C^J (sin p0, 0*) sin qB f ) = -J v {kR)Y v {kR)b pg -J p {kR)J (l {kR) cos [(tt/2) (p-0]<r«(Jto), (5.3) 



^(^=2 S [F J ,_ ff (2nia)-F J , +ff (2wAa)] > (5.4a) 

71-1 

<r#(*»)=2 S (-)"[F,. f (iiJte)-F f+f (iito)]. (5.4b) 

71=1 

We note that the above expressions are needed only for positive integer indices p, q that 
are both even or both odd. We observe also that the form of the expressions is such that 
those with unequal subscript indices can be evaluated simply in terms of those with equal 
subscript indices. Further, in view of the easily-verified relation 

<rg (ha) = 2< (ka) - <r<» (ka/2) , (5.5) 

values or formulas for the cr's "of the second kind" can be derived readily from the same for 
the as "of the first kind." Thus, the evaluation of the G^, which, on the face of it, originally 
involved double integrals with ka, kR, p, g, and h as independent parameters, has been reduced 

essentially to the evaluation of sums of the form ^ZY 2s {nx) with 7r<><C27r or 2t<Cx<C^tt and 

11=1 

s=Qj 1, 2, . . .. It should be observed, however, that these sums converge extremely 
slowly: for the accuracy desired in the present circumstances the number of terms required 
would be of the order of 10 10 or more, at least for the most important o-'s. Needless to say, 
some method of evaluation other than direct summation is indicated. 

5.2. Evaluation of the Bessel-Function Series 

Two methods of evaluating the series defined in (5.4) will be discussed briefly here. 

a. Conversion to Definite Integrals 

An evaluation of the desired series may be obtained in terms of definite integrals as follows. 
We start with the series of Hankel functions (of the second kind) 

Sx(x)=jr;Hx(nx), (5.6) 

77 = 1 

where x is real, positive, and not an integer multiple of 2x. For the Hankel functions we 
use the integral representation [11] 

H x (x)=- f e-* ito M-*tfr 

71" JC 
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where £=£-\-jri and C is a suitable path, to be chosen in a moment, between the limits — j^> 
and 7T+JOO . With the aid of the well-known expression for the sum of a finite number of terms 
of a geometric series, one finds that the M-term partial sum of (5.6) can be written 



M 

J2H x (nx) 






(5.7) 



Among the infinitely many poles of the integrands in these integrals, we note in particular 
£^|and £ TO =x— %' m , where i m is the principal value of sin ~ l (2rriT/x), m = 0, 1,2,.. ., Q, and Q 
is the largest integer such that 2Qt/x<C1. The path C is now chosen to comprise (except for 
indentations) the negative imaginary axis, the real axis between and ir, and the upper half of 
the line ir+jrj; indentations, of radius p, are made at each of the above identified poles in the 
manner shown in figure 7, which illustrates C for the case Q=l. Taken on this path, the 
second integral in (5.7) approaches as M increases indefinitely; as a detailed examination 
shows, this happens because on the various parts of the path the integrand approaches or it 
becomes increasingly oscillatory (or both) as M— ><» . Hence we have 



S: 



T JC ' 






This result can be expressed in terms of real integrations plus contributions from the 
indentations. In the limit as />— >0, the indentations yield 



R '°=h 



T?" 



(5.8a) 



)V 



r * j oo 



V^J 



J^XA 



t: 



-joo 

Figure 7. Path of integration. 



and for m > 1 , 



R'm- 



pM, 



V^ 2 -(2m7r) 2 



Rm = 



e ^ (ir ~U 



^x 2 -(2m7r) 2 



(5.8b) 



Here the notation corresponds to that used for the respective poles and the factor 1/ir standing 
in front of the integral is included. Provided X is an even integer, the straight-line segments of 
the path of integration yield* 



-- <5 oX +lim - 



>-^0 7I\L Jp [ 



cosh Xr? tit) D f ^ 2 . /l . A ,J] 

"^^1 J /> C ° S? \2 X Sin? / f j 



(5.9) 



where 8 \ is the Kronecker delta. It should be noted that when the integrand in the finite 
integral is singular, the Cauchy principal value is to be taken, as indicated by the symbol P 
(the indicated limit as p^O essentially is also a Cauchy value). 
Upon comparing (5.4) and (5.6) it is seen that 

erg (to) = 2 Im [S 2p (2ka)-S (2to)]. 
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It happens that for \=2p (where p = 0, 1,2,. . .) the net imaginary contribution of (5.8) to 
S 2p vanishes. Consequently (5.9) yields 

< (ka) =lj° ^?JT_!_\ +t PJ*'** >' 2 Pi ^t (ka sin Q « (5. 10) 

wherein the integrands are now well-behaved at their lower limits and thus we were permitted 
to put p = 0. The corresponding result for a™ is 

<4W=-|J7 7^ v +i ~t p j7 sin2 ^ ta " (y sin m ' (5 - n) 

which is easily obtained from (5.10) with the aid of (5.5). 

(An interesting by-product of the above development, obtained by taking the real parts of 

(5.6), (5.8), and (5.9), is 

1 ;; i ^ 7 ( n A ] i o^h (>08 t 2 P sin " 1 ( 2m ^)l 
r } o op +ZjJ2 P {rix)—-+22_ l , 

wherep=0, 1,2,. . .; ./• and Q are as specified above. This result was obtained by Ignatowsky 
[12] by a different method.) 

Thus far numerical results for the o-'s have been obtained primarily from (5.10). The 
integrals were evaluated numerically by Paul F. Wacker and William W. Longley, Jr., using 
Legendre-Gaussian quadrature formulas programed for machine computation. 6 Values 
for o { v \l were calculated from values of a™ using (5.5). Results for the two lowest values of 
j) are given in table 2. (More extensive tables are planned for the previously-mentioned 
separate publication.) 

b. Conversion to More Rapidly Convergent Series 

An alternative method of evaluating the o-'s is afforded by a transformation of the Bessel- 
function series involved into more rapidly convergent series. This was accomplished a good 
many years ago by W. von Ignatowsky [12] in connection with his work on diffraction by 
gratings — a subject that will be recognized as being rather closely related to the subject of the 
present paper. Ignatowsky obtained 

* ± r„(w0=log ~+£ 1+ S (l- T J=\ (5.12a) 

71=1 7 n=\n n = Q+l\n -yjn 2 K 2 — 1/ 

f-V (nr) 1 ^ (2/c) 2 '(ff + M-l)!£ 2 M 

n^i^nx)-^ 2 2_ { (2ti\{p-y)\ 

^ sin \2p sin _1 (m)] , . ~^ (nK—-y/n 2 K 2 —l) 2p 
-*2J- (-) p K 22 / 2 2 i (5.12b) 

n=l \l—7l 2 K 2 n=Q+l -vVk — 1 

In these expressions k=2tt/x; #=1,2,3, . . . ; Q is as defined above; log 7=0.57721566 ... is 

the Euler-Mascheroni constant; and the B 2ll are the Bernoulli numbers, here so labeled that 
J3 2 =l/6, J8 4 =l/30, ^=1/42, # 8 =l/30, # 10 =5/66, etc, Equation (5.12a) is fairly well known 

00 

at the present time. (A few numerical values are given for the series ]C (—) n Y (nx) in 

71=1 



6 The estimates of accuracy of computed quantities given in this paper were made by Mr. Longley, who was responsible for most of the detail 
of the calculations. 
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[13].) The series for a[f obtained from (5.12a) and (5.12b) has been duplicated by an in- 
dependent mode of derivation; this represents a partial check of the complicated expression 
(5.13b). Thus far in the present work, Ignatowsky's series have been used only to spot check 
some of the numerical results obtained from the definite integrals. The series are especially 
suitable for either estimation or calculation of the o-'s when p is large. 



The author has benefited from the discussion of many points with Paul F. Wacker and is 
indebted to David F. Wait for assistance with some of the theory and computations. 
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Selected Abstracts 

A comparison of experimental and theoretical rela- 
tions between Young's modulus and the flexural 
and longitudinal resonance frequencies of uniform 
bars, S. Spinner, T. W. Reichard, and W. E. Tefft, 
J. Research NBS 64A, No. 2, 147 (1960). 

The relations between Young's modulus and the mechanical 
resonance flexural and longitudinal frequencies of two sets of 
steel bars has been established by an empirical method. 
The agreement between theoretical and experimental curves 
for flexural vibrations appeared somewhat better than that 
for longitudinal vibrations. 

On the mode theory of very-low-frequency propaga- 
tion in the presence of a transverse magnetic field, 
D. D. Crombie, J. Research NBS 64D, No. 3, 265 

{I960). 

The effect of a purely transverse horizontal magnetic field on 
the propagation of VLF waves is considered. It is shown 
that the magnetic field introduces non-reciprocity, and that 
for propagation along (he magnetic equator, the rate of 
attenuation is less for west-to-east propagation than for east- 
to-west propagation. 

On the theory of reflection of low- and very-low- 
radiofrequency waves from the ionosphere, J. R. 

Johler and L. C. Walters, J. Research NBS 64D, No. 

3, 269 (1960). 

The rigorous application of the magneto-ionic theory to the 
calculation of reflection coefficients for a sharply bounded 
ionosphere model is discussed. This is a comprehensive 
paper on the details which pertain to the rigorous evaluation 
of the reflection coefficients. The paper is illustrated with 
computations applicable to the /)-region or the i^-region of 
the ionosphere. The quasi-longitudinal approximation is 
derived from this theory and the range of validity of this 
approximation is illustrated. The restrictions imposed by 
the use of a sharply bounded model ionosphere are discussed. 

Focusing, defocusing, and refraction in a circularly 
stratified atmosphere, K. Toman, J. Research NBS 

MD, No. 3, 289 (1960). 

Focusing, defocusing, astronomical refraction and path length 
of rays as a function of the departure angle A of the ray at 
the source is described for cases with the source outside, 
inside, or on the boundary of a circular stratification. Rel- 
ative to zero elevation angle symmetrical and centrosym- 
metrical distributions are found. 

Impedance characteristics of a uniform current loop 
having a spherical core, S. Adachi, J. Research 
NBS 64D, No. 3, 297 (1960). 

The radiation impedance is derived by the EMF method in a 
convenient form as the sum of the self-radiation impedance 
of a loop in the free space and an additional term due to the 
reaction between the loop and the sphere which is propor- 
tional to the well-known expansion coefficient of a magnetic- 
type scattered wave from a sphere in an incident plane wave. 
The first antiresonance frequency has been given in the form 
of a universal curve for a very small uniform current loop 



with core of an arbitrary composition of ju. s and e s , subject 
to the condition that the refraction coefficient N=sJfjL S e s is 
extremely large. Some numerical calculations show that 
high-^t core is desirable for a comparatively lower frequency 
region, and high-e core is rather desirable in an antiresonance 
region. 

Basic theorems in matrix theory, M. Marcus, NBS 

Applied Math. Series 57 (1960) 15 cents. 

This is a survey of the basic identities and inequalities of 
matrix theory. Included are results dealing with elementary 
properties, canonical forms, invariance, congruence, com- 
mutativity, orthogonalization, eigenvalues, determinants, 
submatrices, rank, determinant and rank inequalities, nu- 
merical methods for inversion and eigenvalues, condition 
numbers. 

Relaxation processes in multistate systems, K. E. 
Shuler, Phys. of Fluids % No. 4, 442 (1959). 

An analysis is made of the relaxation of "multistate" systems, 
i.e., systems with more than two quantum states or two 
different chemical species, for linearized processes described 
by the " Master Equation" of the theory of transport processes. 
These results are compared with those obtained from the 
analysis of two-state relaxation processes and the concept of 
"relaxation time" is discussed in this framework. A dis- 
cussion is presented of the transformation of microscopic to 
macroscopic relaxation equations. The existence of periodic 
and/or aperiodic oscillatory solution of the linear multistate 
relaxation equation is investigated. It is shown that the 
multistate relaxation equations admit of aperiodic oscillatory 
solutions. 

Invariant and complete stress functions for general 
continua, C. Truesdell, Arch. Rational Meek, Anal. 

l 9 No.4, 1 (1959). 

This paper presents an organized survey of results, with 
proofs, concerning exact general solution of the underdeter- 
mined equations of motion of a continuous medium. Special 
attention is given to a flat space of n dimensions, to an 
n-dimensional space of constant curvature, and to a 2-dimen- 
sional surface. Two methods are used: (1) the representation 
of a solenoidal vector as the curl of a vector potential, and 
(2) a variational formulation, essentially the converse of the 
principle of virtual work. 

The construction of Hadamard matrices, E. C. Dado 
and K. Goldberg, Mich. Math. J. 6, 247 (1959). 

It is proved that a Hadamard matrix of order An can be 
constructed if there exists a transitive permutation group of 
odd order and degree 4n-l whose subgroups leaving one ele- 
ment fixed have three transitivity sets each. 

Bending and stretching of corrugated diaphragms, 

R. F. Dressier, ASME Trans. 81D, No. 4, 651 (1959). 

Solutions of the exact linear elastic shell equations for all 
stresses and displacements in a typical corrugated diaphragm 
are presented for three cases over a 9 to 1 range of thicknesses. 
Results were obtained by numerical integration in an elec- 
tronic digital computer. The effect of thickness variation is 
discussed with respect to both stresses and resultants, and 
peak values needed for design purposes are presented. Cir- 
cumferential and meridional stresses are found to be equally 
important throughout the thickness range analyzed. Bending 
and membrane stresses are likewise equally important 
throughout the range. Peak values in some cases occur near 
the outer rim. 
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Modular forms whose coefficients possess multi- 
plicative properties, M. Newman, Ann. Math. 70, 

478 {1959). 

Let q be a prime, r and s nonzero integers, and 77 (r) the 
Dodekind modular form. All modular forms of type 
97 r (r)r; 8 (gr) such that the related Dirichlet series possesses 
an Euler product are determined. 

The analysis of latin squares with a certain type of 
row-column interaction, J. Mandel, Technometrics 1, 

No. 4, 379 (1959). 

A serious limitation in the use of latin squares is the con- 
founding of the main effect of each factor with the interaction 
of the remaining two factors. In some cases, the interaction 
of rows and columns can be expressed as a multiplicative 
term of assigned factors associated with rows and columns. 
The analysis of such designs is presented in detail. An 
example is discussed. The method is extended to the case 
where the treatment effects are linearly related to a given 
concomitant variable. A brief discussion is given of the 
relation of the proposed method with other tests for non- 
additivity. 

A note on algebras, A. J. Goldman, Am. Math. Mo. 
66, 795 (1959). 

The usual sufficient conditions that both regular representa- 
tions of an algebra A be faithful are shown not to be necessary. 
A necessary and sufficient condition (viz, that no nonzero 
member of A be either a left annihilator or a right annihilator) 
is derived. 

Linear differential equations of the second order with 
a large parameter, F. W. J. Olver, J. Soc. Indust. 
Appl. Math. 7, 306 (1959). 

An investigation is made of the differential equation 

g?=Mp(Z)+ a (Z)}« 

for large values of the parameter u. 
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